arXiv: 1506.07301vl [nlin.CD] 24 Jun 2015 


Averaging Theory for Weakly Non-linear Oscillators 

Aritra SinhaQ 

Department of Physics and Astronomy,National Institute of Technology Rourkela Address: Sector - 2 Rourkela, Odisha-769008 

(Dated: 15 July 2014) 

Here I have first discussed how averaging theory can be an effective tool in solving weakly non¬ 
linear oscillators. Then I have applied the technique for a Van der Pol oscillator and extended the 
stability criterion of a Van der Pol oscillator for any integer n(odd or even). 


WEAKLY NONLINEAR OSCILLATORS 

We consider systems of the form 
X -b X -b eh(x, x) = 0 

where 0 < e <C 1 and h{x,x) is an arbitrary smooth 
function. 

Since such equations represent small perturbations 
of the linear oscillator x + x = Q , they are called weakly 
nonlinear oscillators. Some well known examples are 
the Van der Pol equation 

X + X + e{x^ — l)i; = 0 

and the Duffing equation 

X -\- x + ex^ = 0. 

For the Van der Pol equation in (x, x) phase space, if we 
choose an initial condition close to the origin and e near 
0, the trajectory is a slowly winding spiral; it takes many 
cycles for the amplitude to grow substantially.Eventually 
the trajectory asymptotes to an approximately circular 
limit cycle whose radius is close to 2. 



There are various ways one can exploit the fact, that os¬ 
cillator is ’close’ to a simple harmonic oscillator, to pro¬ 
duce useful approximations. Here I will use the method 
of averaging theory to predict the period and radius of 
the limit cycle. 


AVERAGING THEORY 

Introduction 

Our system is 

x-\-X + eh(x,x) = 0 where 0 < e <C 1 
Let 

x = y 


y = -X- eh{x,y) 

When e = 0, the solutions are 

x{t) = r cos {t + (f ) 

y{t) = —r sin {t + (f) 

where r,^ are constant for simple harmonic oscillator on 
trajectories. 

When e 0, we expect a very slow drift of r and </>. 
They would evolve as expected; might approach a limit 
cycle.Our task at hand is to find the evolution equations 
of the amplitude and the phase given the effect of this 
nonlinear term h{x,y) .We can be definite that since e is 
extremely small, the trajectories would be nearly circular 
and they would have periods of 27r approximately. 

Now let 

x{t) = r(t) cos (t -b 
y{t) = -r{t) sin {t + cj){t)) 

Note here that although it might seem that the calculus 
done is wrong in computing y, but it is not so. Instead 
let it be a definition. 

So view this as definition of r{t) and 4>(t) i-e. 

r{t) = \/a:2(t) +y'^{t) 
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Calculating r and (j) 


Main Concept of Averaging 


Let us find equations for : 

2 2,2 
r = X +y 

rf = XX + yy 

= x{y) + y{-x - eh) 

= -e-yh 

= —eh{—r sin {t + (j))) 


r — eh sin (t + 4>) 


d , ,, ,, d f y(t) \ 


1 + — 


x{t)J 

d{-l) 


1 + S / dt 


x^ + y^ 

+ y"^ 

1 

x^ + y"^ 

1 


, 1 . y ■ 

-\ -y — 2 ^ 

X x^ 


y ■ y 

X^ X 


[yx - xy] 


= 1 


x^ + y"^ 
ehx 


+ exh\ 


x^ + y^ 


Now since 


Now what we have to do is exploit the separation of 
time scales of the evolution equations - fast oscillation 
versus slow drift. The main idea is to iron out the fast 
oscillations by averaging over one cycle of length 27r (1 
oscillation) so that we can find out the explicit depen¬ 
dence of the evolution equations on e. 

A practical example of the averaging theory would be 
our primary detector, eyes. All the things we see around 
us are made up of rapidly vibrating atoms. But our de¬ 
tectors average out this fast oscillation and we see only 
the slow scale dynamics explicitly. 

Given g{t),\et us define the average over one cycle about 
the point t as 

2 ^ rt+TT 

< 5 > W = ^ J 9 is)ds. 

So let us write out the equations for f and ej). 
r =< eh sin {t + (j)) >t 


(j) = (^ — cos {t + (j)) 

These equations are exact since R.H.S. depends on 
not f, (j). 

Now here you will notice a subtle difficulty, since we want 
a dynamical system with f and 0 so that we can com¬ 
pute the averages.We can get over this problem with an 
approximation.We well know that there is only a slight 
difference between r and f. [It is comparable to our height 
measured right at this instant versus our height averaged 
over the past. That is sort of the same thing with a slight 
error.] It is here that we are going to introduce an ap¬ 
proximation. 

Over one cycle, 


r = r + 0(e) 


x^ +y‘^ = and x = r cos {t + 

we have. 


* 

(p — — cos {t + (p) 
r 


So we get slowly varying amplitude and phase, as ex¬ 
pected , of the order e. 


0 = 0-1- 0(e) 

Now we will replace r,0 by f,0 in the evolution equations 
which causes O(e^) in the ODE’S. The beauty of this is 
that we get autonomous equations which we can analyse 
easily with phase plane methods. So the formal evolution 
equations are 

f —< eh sin {t + <P) > +0(e^) 


f = 0(e), 0 = 0(e) 




cos (t + <p) 


+ 0(e2) 
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We should treat r and (j) as constants while performing 
the averages. 


= —er” (cos”''"^ 9 sin 6) + e (sin 9 cos 9) + O(e^) 

= 0(e2) 


APPLICATION ON GENERAL VAN DER POL [•.■ both (cos”+^ 6>sin6») and (sin6»cos6>) are equal to 0.] 

OSCILLATOR 


Now we shall apply this averaging technique to a gen¬ 
eral case of Van der Pol’s oscillator.The general Van der 
Pol equation is- 

x-\-x-\-ex(x^ — \) =0, 0 < e <C 1 'n' is a natural number 

I will try to show that when n is odd, r is going 
to grow exponentially and when n is even, r approaches 
a limit cycle. 

I am going to use ”< >” symbols which essentially 

means average over one cycle. We need to remember 
from previous discussions that x = r cos (t + (f),]! = 

—r sin (t -b 4>) and x = y. Let {t + (f) be replaced by 
9 in the following math. [Now let 


When n is Odd 

Let n = 2m -b 1. 

(cos^^+i^sin^d) 

9 9 d9 

t = sin 9, then dt = cos 9d9] 



h{x, x) = (x” — l)x = (x" — l)y 


= / (1-P 


dt 


= (r" cos" (t -b 0) — 1)(—r sin {t + <p)) + 0(e) 

= —cos" 0 sin 0 -b f sin0 -b 0(e) 

f =< eh sin 9 > -bO(e^) 

= (e(f sin0 — cos" 9 sin 9) sin0) -b O(e^) 

= ef (sin^ 9) - er"+i (cos" 9 sin^ 9) -b ©(e^) 

= Y (cos" 9 sin^ 9) + ©(e^) 

0 = ^ — cos (t -b 0)^ -b O(e^) 

= cos" dsin0 -b f sin0) cos0^ -b O(e^) 

= (e(—f" cos""*"^ 0 sin 0 -b sin 9 cos 9)) + O(e^) 


= 0 . 

So for odd n, 



If we neglect O(e^) error, then 

df er 
dt 2 



=> logr=-t-bA: [where k is any constant] 

. r = pexp (^-t^ [where p is a constant of integration]. 

So the final approximate solution for odd n is 
X = pexp cos (t -b 0). 

Therefore we assert that x is going to blow up given suf¬ 
ficiently long time and the phase is going to change neg¬ 
ligibly. 
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When n is Even 

We can be sure, that since cos" 9 sin^ 0 is > 0 through¬ 
out, (cos" 9 sin^ 9) must be equal to some positive quan¬ 
tity say k"^. 


r = - er k 

2 

= f(l-2r-fc2) 

So f has one real positive root Whenever f is 

greater than ( 5 ^) " ,~r is less than 0 and f continues to 
shrink till it reaches {^) " • And whenever f is less than 
( 2 S 2 )"; ^ is greater than 0 and f continues to grow till 

it reaches ( 2 ^ 2 )"- So this means that the oscillator is 
going to approach a limit cycle of amplitude 2. This is 
the beauty of the problem when n is even. Now let’s do 
a problem based on this. 

Example: Let us consider a Van der Pol oscilla¬ 

tor 


Plot of time-derivative of r with r 



X -I- a: -I- ex(x^ — 1) = 0, 0<e<Cl 


h{x, x) = (x^ — l)x = (x^ — l)y 


= (r'^ cos^ {t + 4') — 1 )(—r sin (t -|- (j))) + 0 (e) 


r =< e/isin (t + (j)) > -|- 0 (e^) 


= (e(r sin 9 — cos" 9 sin 9) sin 9) + O(e^) 
= ef (sin^ 9) - ef^ (cos" 6 lsin^ 9) + 0{e^) 


er er 




+ 0(e2) 


er 


= V(4-r")+0(e") 


From the diagram it is clear that as t tends to 00 , 
value of f tends to 2. It means that after sufficiently long 
time,the Van der Pol Oscillator settles down to a limit 
cycle of amplitude 2 . 

(j) ^ O(e^) [proved earlier in the general case] So ej) 
changes on super slow time-scale. Period of the Van der 
Pol oscillator = 27r -|- O(e^). 
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